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Abstract 

A moving boundary separating two otherwise homogeneous regions of a dielectric is known to 
emit radiation from the quantum vacuum. An analytical framework based on the Hopfield model, 
describing a moving refractive index step in 1+1 dimensions for realistic dispersive media has been 
developed in [1], We expand the use of this model to calculate explicitly spectra of all modes of 
positive and negative norm. Furthermore, for lower step heights we obtain a novel set of mode 
configurations. This leads to a realistic emission spectrum featuring black hole- and white hole 
emission for different frequencies. We also present spectra as measured in the laboratory frame 
that include all modes, in particular a dominant negative norm mode, which is the partner mode in 
any Hawking-type emission. We find that the emission spectrum is highly structured into intervals 
of emission with black hole-, white hole-, and no horizons. Finally we estimate the number of 
photons emitted as a function of the step height and find a power law of 2.5 for low step heights. 
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I. INTRODUCTION 


In an effort to tie the theory of quantum helds together with that of general relativ¬ 
ity, Hawking showed [2, 3] that black holes emit a steady thermal flux: in contradiction 
with their name, black holes are not black. The phenomenon, now called Hawking Radia¬ 
tion (HR), arises from quantum vacuum fluctuations at the event horizon. It connects the 
realms of quantum physics, Einstein’s theory of relativity, and thermodynamics. However, 
it is essentially impossible to observe in astrophysics because of its low temperature, which 
depends on the surface gravity at the horizon, i.e. the mass of the black hole. A solar-mass 
black hole has a temperature well below that of the cosmic microwave background. 

Fortunately, Unruh [4] showed how wave media, if moving, can be used to realize gravity 
analogues. A flowing fluid exhibiting a gradient of flow velocity from subsonic to supersonic 
can mimic the flow of space at the event horizon: there is a boundary, beyond which an 
acoustic wave propagating against the flow will be carried downstream. This is the analogue 
of an event horizon. In complete analogy, quantum hydrodynamical fluctuations in a moving 
fluid are thus predicted to convert into pairs of phonons at the sonic horizon. Analogue 
systems produce analogue HR, providing a realistic chance to observe this particle emission 
from the quantum vacuum [5, 6]. 

Until recently, proposed analogue systems to test Hawking’s prediction by experiments 
suffered from a low characteristic temperature [4, 7-13]. In 2008, Philbin et al [14] demon¬ 
strated the feasibility of creating artificial event horizons with a moving refractive index front 
(RIF) in dispersive optical media. They estimated a temperature of a thousand Kelvins and 
ushered-in the held of optical analogues [1, 15-17]. The idea behind optical event horizons 
is to create a change in the refractive index, i.e. the speed of light, with light itself. For 
example, a pulse of light modihes the index by the optical Kerr effect. Light under the pulse 
will be slowed and thus the front of the pulse exhibits -for some frequencies- a black-hole 
type horizon capturing light. The back of the pulse acts as an impenetrable barrier, a white- 
hole horizon. The event horizon separates two discrete regions: under the pulse, where light 
is slow and the pulse moves superluminal, and outside the pulse, where the pulse speed is 
subluminal. 

To date, optical HR is yet to be discovered. Detailed analytical predictions of the wave¬ 
length and intensity of the radiation for an actual experiment are difficult to make due to 
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the extreme conditions at the horizon. In the above example, the length of the pulse front 
or back has to be comparable to the wavelength of radiation [14]. Nonetheless, progress has 
been made towards an understanding of the critical conditions needed to observe optical HR 
[1, 18-25]. In particular, Finazzi and Carusotto [1, 18, 19, 25] introduce a fully quantized 
analytical model for spontaneous emission at a moving boundary in a nonlinear dielectric. 
At this boundary between two multi-branch dispersive media certain modes may experience 
either analogue black- or white hole or horizonless conhgurations, leading to mode mixing 
and spontaneous emission of radiation. All conhgurations lead to the spontaneous emission 
of radiation due to the mixing of modes with opposite norm. In order to maintain these 
different conhgurations, they hnely adapt the velocity of the RIF when changing the nonlin¬ 
earity. They discover that emission is dominant over optical frequencies. Hence they focus 
on emission spectra in positive norm optical modes only. All of quantum vacuum emission is 
accompanied by emission of negative norm modes, which play the role of the partner mode 
in the Hawking ehect. This is relevant, in particular, because these modes emit at diherent 
laboratory frequencies compared to their positive norm partner modes. Furthermore, the 
nonlinearity in the experiment typically changes independently of the RIF velocity, lead¬ 
ing to a very diherent spectral structure as well as scaling of the signal with nonlinearity/ 
intensity. 

In the present paper we use the model [1] to reveal these properties of the quantum 
vacuum emission. We hrst expand the model to consider emission from all positive and neg¬ 
ative norm modes at any frequency and change of refractive index. Hence we obtain emission 
spectra for diherent index changes without tuning of the RIF velocity. This corresponds to 
most experimental situations where the velocity of a RIF would not change with height. 
A resulting consequence is that the mode conhgurations can be more complicated than for 
a black- or a white hole. In fact, we obtain spectra which combine black- and white hole 
emission. The system is then an object that simultaneously radiates as analogue black and 
white hole. We convert the spectra to the laboratory frame, including all mode contribu¬ 
tions, inclusive of the important negative norm ones. The lab spectrum is characteristically 
structured into intervals of emission with black hole-, white hole-, and no horizons. Finally, 
we observe that the total photon hux associated with black hole emission rises with a power 
of about 2.5 for low refractive index changes, whereas the emission bandwidth increases 
linearly. For larger index changes, Sn > 0.05, which are difficult to reach experimentally. 
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the rise slows to a lower power law as the number of participating modes saturates. These 
novel features will be useful to identify emission in future experiments on optical Hawking 
radiation. 

We hrst introduce the theoretical model of the scattering of vacuum modes at the horizon. 
We explain how the interaction of light and matter in a uniform dispersive medium is 
modeled. We identify the eigenmodes and study their properties and then introduce the 
moving Refractive Index Front (RIF) separating sub and superluminal regions. We construct 
eigenmodes of the non-uniform medium with the RIF and describe the scattering/ mode 
conversion process at the RIF by the scattering matrix. We then quantize the field modes 
and calculate the photon flux density in the moving and laboratory frame. Finally we 
compute the spectra of emission from bulk fused silica in both frames. We also integrate the 
spectra to evaluate the total emission and its scaling with the refractive index height. The 
spectra allow us to indentify in detail the contributions of the various modes to the emission 
and the role of analogue event horizons. 


II. LIGHT MATTER INTERACTIONS IN A DISPERSIVE MEDIUM 

We start by laying out the canonical quantization scheme for the electromagnetic held in 
the linear and non-uniform dielectric. We restrict ourselves to a one-dimensional geometry 
and scalar electromagnetic helds and operate at frequencies sufficiently far from the medium 
resonances to neglect absorption. Our analysis is based on the Hopheld model [26] and 
we follow the treatment presented in [1] before expanding it to include all modes at any 
frequency and refractive index change. 

The step in refractive index (RIF) is propagating in the positive x-direction and the 
electric held is orthogonal in the z-direction. Light propagation in the homogeneous medium 
is sufficiently well described by the interaction with a set of three polarization helds, which 
phenomenologically reproduce the dispersion relation of transparent dielectrics. We state 
the Lagrangian density C describing the coupling of the vector potential A{x,t) with the 
three polarization helds Pi, harmonic oscillators of elastic constant and inertia ^ , 

in the laboratory reference frame: 

{djAf {d^Af ^ f (dj. 

Sttc^ Stt \ 2 Ki 
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The term linear in A in Eq. (1) describes the coupling between the helds. This Lagrangian 
density accounts for the free space and medium contributions to the held through the hrst 
two terms and the sum, respectively. Dispersion enters as a time dependence of the addends 
of the summation. 

It is convenient to express the Lagrangian density in the moving frame by applying a 
Lorentz boost. The system is stationary in the reference frame comoving with the RIF at 
velocity u. We denote the space and time coordinates in the moving frame by, respectively, 
( = 'j {x — ut) and r = 'y (t — x^') with the Lorentz factor 7 = (1 — v?/(A) . We obtain 

the Hamiltonian density by varying the Lagrangian density with respect to the canonical 
momentum densities H^i = and Hp. = of the light and polarization helds. From 

there follow the Hamilton-Jacobi equations, the equations of motion for the helds. 

The linear Lagrangian (1) leads to the generic Sellmeier dispersion relation of bulk trans¬ 
parent dielectrics: 




1 + E 


Attk^ 


TT 1- 
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( 2 ) 


where u and k are the laboratory frequency and wave number, respectively. The dispersion 
relation ( 2 ) is plotted in hg. 1 : there are eight branches, four with positive laboratory 
frequency and their four negative frequency counterparts - symmetric about the /c-axis. To 
hnd orthonormal plane-wave solutions, we dehne a scalar product on the set of the solutions 
of our Hamilton-Jacobi equations generalized to complex values. 


{VuV2)='- 

where V is the eight-dimensional vector 


r)vV2{C,r), 


( 3 ) 


V = {A p, P 2 P 3 Hp, Hp, np 3 )^, 


( 4 ) 


and 7] is the selection matrix rj = ( o') with being the 4x4 identity matrix. This 
scalar product is a conserved quantity in r and therefore the total norm of the quantum 
states is conserved. 

It can also be shown that modes with positive laboratory frequency have positive norm, 
and those with negative frequency have negative norm [1, 14, 25]. In other words, modes 
belonging to the upper (lower)-half plane in hg.l have positive (negative) norm. In the 
comoving frame, positive frequency waves with negative norm appear, which are associated 
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with spontaneous emission from the quantum vacuum. These generated waves have negative 
laboratory frequency. Negative norm modes were recently observed in water experiments 
[8, 27] and in optics [22, 28]. Due to the conservation of norm, the generation of negative 
norm waves signihes a simultaneous increase in positive norm waves, the generation of 
correlated waves. 


III. MODE CONFIGURATIONS AT A REFRACTIVE INDEX FRONT 


In the previous section we presented a canonical model of light and matter interaction 
in a dielectric medium. In particular, we highlighted the existence of positive and negative 
frequency modes. We now extend this model to a non-uniform medium. 

The simple geometry of a refractive index front (RIF) is shown in £g. 2 in the comoving 
frame. The medium is composed of two homogeneous regions, separated by the RIF at C = 0 
creating a step in the refractive index. The boundary at C = 0 constitutes an infinitely steep 
RIF which propagates in a steady and rigid way in the positive ( direction. The right (left) 
region has dispersion parameters KiR^HiL) and XiR^Xn)- Experimentally, this front could be 
realized by an asymmetric intense pulse of light with a very steep front on one side, similar to 
a self-steepened soliton pulse [29]. The front is approximately infinitely steep if it is shorter 
than a wavelength. The step height is the refractive index change between the regions, Sn, 
defined by 

n (C) =nL0 (-C) +nR9 {() = nR + Sn9 {-() (5) 


and illustrated in fig. 2; 9 {() is the Heaviside step function. The index change is described 
by the scaled Sellmeier coefficients PiR = cr/dij? and = crX'fj^, where: 

2nR Sn 


(j ~ 1 -|- 


nj- 1 


( 6 ) 


and Ur is the refractive index on the right side [1]. 

In order to identify the allowed modes of propagation, we consider plane wave solutions 
subject to the dispersion relation in both regions for a given comoving frame frequency 
u' = 'y{u — uk). This is useful, because the frequency u' is conserved, an expression of 
energy conservation in the moving frame [1, 14, 25]. Solutions of fixed u' are found at the 
intersection points between a line of constant u' with the various polariton branches in the 
dispersion diagram (red circles in fig. 1). Note that we only consider low positive comoving 
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frequencies u' without solutions in the uppermost branch. At any frequency u' there exists a 
set of 8 discrete (ca, k) solutions to the dispersion relation (2). We either hnd 8 propagating 
modes or 6 propagating modes and two exponentially growing/decaying modes, that take 
on complex u and k. 

Emission spectra with 8 propagating modes on only one side of the boundary were cal¬ 
culated in [1, 19]. Here we address the experimentally relevant case allowing for 8 modes 
to propagate on either side of small refractive index changes. In particular, consider the 
positive frequency optical branch (where A 2 < A < A 3 ) in the moving frame, depicted in Eg. 
3. The black (orange) curve is the branch on the right (left) side of the RIF. Again, the num¬ 
ber of solutions depends on u'. On either side of the RIF there is at least one propagating 
optical mode for all u'. There also is a frequency interval where three propagating modes 
exist: between the two horizontal dashed black and the two dashed orange lines in Fig. 3, 
corresponding to the right and left region, respectively. Hereafter these frequency intervals 
on either side of the RIF are referred to as the subluminal intervals (SLI) and 

Inside a SLI, one of the three mode solutions has a positive comoving group velocity 
This unique mode allows light on the right of the RIF to escape from it. This middle 
optical mode (cf. hg. 1) on the right is called moR in what follows. The other two 
modes have negative comoving group velocity; they move into the boundary from the right. 
There is a lower (upper) optical mode loR {uoR). On either side we can order the modes 
by the comoving wave number k' and obtain (cf. hg. 3). In 

the laboratory frame this translates into (hg. 1 ). All three 

optical-branch modes of the SLI have positive laboratory group velocity on both sides of the 
boundary: to an observer in the laboratory frame they propagate in the same direction as 
the RIF. 

Beyond this interval (i.e. u' ^ [i^min)'^max]) propagating mode remains. Two 

complex-wavenumber roots emerge as pairs of exponentially growing and decaying modes 
that do not propagate. For u' < only mode uoR/L remains a propagating mode whereas 
for oj' > only loR/L remains. It can also be seen in hg. 3 that for all u' there is one 

propagating mode that belongs to the negative optical-frequencies branch. This mode has 
negative norm and will hereafter be referred to as mode noR/L. 

The subluminal frequency interval of the left region is in general diherent from that of 
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the right region. Therefore, there exist 5 different combinations of modes across the RIF, 
also shown in £g. 4: 

1. ca' < one optical propagating mode {uoR/L) exists (real u and k) on either 

side of the boundary. All propagating modes exhibit negative group velocities in the 
moving frame and therefore no optical horizons exist. 

2. < u' < three propagating modes {loL, moL, and uoL) on the left of the 

boundary exist, but only one mode {uoR) on the right. Only mode moL has positive 
group velocity on the left. Light in this mode experiences a white hole horizon at the 
RIF as it can approach, but not enter the right region. To our knowledge, this mode 
configuration has not yet been described. 

3. < u' < three propagating modes {loR/L, moR/L, and uoR/L) exist on 

either side of the boundary. Modes with negative and positive group velocity exist on 
either side of the RIF, therefore the RIF is not a one-way door and no horizon exists 
for waves of this moving frame frequency. 

4. < oj' < WmaxK- propagating mode {loL) exists on the left of the 

boundary, but three modes {loR, moR, and uoR) do on the right. Only mode moR 
has positive group velocity on the right side. Light experiences a black hole horizon at 
the RIF as it cannot propagate to the right from beyond the RIF. The region on the 
left (right) of the boundary corresponds to the inner (outer) region of the analogue 
horizon. 

5- one propagating mode (loR/L) exists on either side of the boundary. All 

propagating modes exhibit negative group velocities and therefore no optical horizons 
exist. 

Note that for very high refractive indices on the left there is only mode loL and no subluminal 
interval in this region and out of the general case only configurations 1 (with loL instead of 
uoL), 4, and 5 remain. In this case, modes within the SLI on the right of the boundary always 
follow configuration 4. The 5 configurations introduced here describe a typical experimental 
combination of horizons with low refractive index changes. 



IV. SCATTERING AT THE RIF 


Having derived solutions on either side of the RIF, we are now constructing ‘global’ 
solutions, i.e. solutions to the equation of motion that are valid in both regions. These 
modes correspond to waves scattering at the RIF and they describe the conversion of an 
incoming held, even if in the quantum vacuum state, to scattered helds in both regions. We 
follow the approach [1, 26, 30] to construct these modes and their scattering matrix and then 
to quantise the solutions to hnd the photon huxes due to spontaneous particle creation. 

The solutions to the dispersion relation correspond to local modes (LMs) that are dehned 
in either the left or right region. We construct ‘global modes’ (GMs) V as: 

V = L‘V^ 0(-() + R'-Vg 0(0, (7) 

a a 

where L°‘ (i?“) describes the strength of mode a on the left (right) side of the RIF. It can be 
shown that across the boundary, vector potential, polarisation helds, and their derivatives are 
continuous [1], which constrains half of the coefficients in Eq. (7). We follow the procedure 
described in [31] and consider GMs whose asymptotic decomposition comprises only a single 
LM with group velocity towards {in) or away from {out) the RIF. Thus there are as many 
of these GMs as there are propagating local modes. Half of the GMs emerge from a dehning 
LM a that moves towards the RIF, forming ‘global in modes’ V““. The other GMs are 
‘global out modes’ if a is a LM now moving away from the RIF. The sets of and 

yout ] 2 iodes are two basis sets of modes, because the LMs are the complete physical (i.e. 
non-divergent) solutions in the asymptotic regions. Hence the scattering matrix S is the 
transformation of modes from the out-basis to the in-basis: 

yin,a ^ ^ ^a/3yout,/3_ 

0 

Scattering and spontaneous photon creation occurs as the input state vacuum does not 
correspond to the vacuum state in the out-basis, i.e. the spontaneous emission follows from 

We postulate the equivalent of the standard equal-time commutation relations on A and 
Pi and thus quantize the local held modes and their momenta: 

WO.nxK'il^wiK-C'). [P.(C),np,(c')] =«'5.,'5(C-C')- (9) 
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We expand the real field V (Eq. (4)) on a basis of local frequency eigenmodes 


V 


duj 






+ d 


af 


( 10 ) 


that are properly normalized with respect to the scalar product (3) under the condition 



5(6^2 ^l) da2ai- 


( 11 ) 


The operators a"/ and are the annihilation and creation operators of the field mode 
a. Alternatively, we can expand the field over positive frequencies only, including negative 
norm modes in the expansion: 


V 




+ ff.c., 


( 12 ) 


where P (N) is the mode set of positive (negative) norms. Writing the global field V in the 
basis of global-in-modes induces quantisation onto the GMs: 


V 




+ P.C. 


(13) 


and a similar expression can be written for V using the global-out-modes. The expansion (13) 
for in- and out modes defines the annihilation and creation operators for the global modes as 
well as the transformation between in- and out creation and annihilation operators. If is 
a row vector containing all the annihilation (or creation) operators for positive (or negative) 
norm global-in-modes and correspondingly for the out-modes, then the transformation 
of operators is [1]: 

^out _ ^in 


Using Eq. (14) we obtain the flux density of photons in mode a, the number of particles 
per unit time At and bandwidth in the moving frame, as [1]: 


j/a _ 



(15) 


In terms of the scattering matrix, the flux density in mode a follows from Eq. (14) and Eq. 



where a is a mode with a norm opposite in sign to a. Spectra for positive norm modes and 
for single SLIs were calculated with this relation in [1, 19]. We now move on to obtain the 
spectra of emission for any frequency in all modes, all mode conhgurations (see sec. Ill), 
and for a variety of refractive index changes Sn. 


V. EMISSION SPECTRA AND PHOTON FLUX 

We proceed to use the S matrix to compute the spectra of emission into all optical modes 
as seen from the moving and the laboratory frame. We consider modes in bulk fused silica. 
The material resonances are Ai_ 2,3 = 9,904nm, 116 nm and 68.5 nm, respectively, and the 
elastic constants are ^ 1 , 2,3 = 0.07142, 0.03246, and 0.05540, respectively. The velocity of the 
RIF is M = 0.66 c, corresponding to a group index of 1.5. 

We hrst consider spectra in the moving frame. Fig. 5 displays the spectrum of emission 
into moR, the unique rightgoing mode (black hole emission), for different index changes 5n. 
Spectral emission is constrained to the right SLI '^maxi?]; where the mode moR exists. 

An optical horizon, however, does only exist for part of this interval, i.e. '^maxi?]) 

because at lower frequencies the SLIs of the left and right overlap (cf. £g. 3). As can 
be clearly seen, the absence of a horizon leads to a significant decrease in the emission, 
i.e. mode coupling, although some emission remains. In case of at a large index change 
(e.g. 5n = 0.056), as calculated and observed in [1], the spectral shape is not thermal as is 
expected from a finite subluminal interval. Here we observe a remarkable feature: the shape 
of the spectrum above where a horizon exists, is independent of the index change. 

This can be more clearly seen in £g. 6 , in which the spectra for lower index changes are 
scaled up to compensate the lower single frequency rate. All traces over orders of magnitude 
of index changes line up to the same shape, making it a universal signature of analogue 
black hole emission. Note also that the shape differs for emission outside the SLI. 

Next we consider spectra of emission into each mode (i.e. noL, loL, uoL, moR) for a fixed 
Sn = 0 . 02 , as displayed in fig. 7. The strongest emission occurs into the optical mode with 
negative norm, noL. This emission is due to coupling with all the other positive norm modes 
in the medium and is strongest where this superluminal mode couples to a subluminal mode 

(^maxL <UJ' < 

It is interesting to calculate the total photon flux over the SLI by inte- 
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grating over the spectrum of fig. 5. In order to convert the flux to a realistic, though very 
approximate photon number, we assume that the RIF propagates over a distance of 1 mm. 
The resulting photon number as a function of index change 5n is given in £g. 9a). The num¬ 
ber of photons excited from the vacuum first grows with power 2.5 of 5n until 5n = 0.052. 
The emission spectrum becomes wider in a linear way, as shown in £g. 9b). Thus the 
emission rate for a single mode increases with This result is surprising and will be 

further investigated elsewhere. For very high induced index change, i.e. for 6n > 0.052, the 
spectral width saturates and the emission rate grows slower accordingly. These index steps 
are difficult to reach experimentally by nonlinear optical pulses though. 

The spectra calculated above would be observable in the moving frame. In an actual 
experiment the detectors are located in the laboratory frame and thus the spectrum observed 
is different. For each mode, the rate of photon production per unit time and unit frequency 
u in the laboratory frame is derived from the moving frame rate as [1]: 

V 

where Vg{u) is the laboratory group velocity at u. The spectral density as a function of 
frequency converts to wavelength by the factor (27rc). A lab spectrum for a single mode 
was calculated in [1]. However, this spectrum is not observable for the following reason: on 
either side of the RIF, each moving frame frequency u' corresponds to up to eight different 
laboratory frequencies depending on the mode involved, as in fig. 1. Note that we calculate 
the emission with positive laboratory group velocity only. Thus each mode has a unique 
direction of group velocity in both frames independent of the side of the RIF. Hence, of the 
modes left and right of the RIF, half are defining global out (emission) modes. Conversely, 
the emission at a hxed laboratory frequency may arise from several optical modes. We 
proceed to calculate for successive frequencies u the emission rates in each associated global 
out mode, obtaining a laboratory frame spectrum for each mode. We then add up the 
spectra to obtain the total spectral emission. Starting at Um, the lowest positive optical 
frequency, we obtain emission into outgoing mode loL on the left, then moR on the right 
followed by uoL on the left and hnally noL on the left. 

Fig. 8 shows laboratory spectra for three index changes. As the spectrum is composed of 
contributions from different modes for different mode configurations, it exhibits a number of 
sharp features. The largest spectral density is obtained around 250 nm and corresponds to 
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emission from the negative norm mode noL. Emission is generated by the pairwise coupling 
of two modes of opposite norm. Hence mode noL is the only negative norm mode on the 
optical branch and covers a rather small laboratory spectral interval, between the violet and 
red lines in fig. 8. Therefore, all emission due to the coupling of two optical modes leaves 
a contribution within this emission peak in the UV spectral range. The coupled positive 
norm mode (i.e. the partner photon), if optical, can be found at the remaining optical 
frequencies. We choose to limit the range of optical wavelengths such that no modes in the 
uppermost branch are excited, resulting in a cut off at 230 nm. Not all coupled mode pairs 
are separated by a black or white hole-type horizon. For example, intervals with horizons, 
as schematically sketched in Eg. 4, are found between the black and orange dashed lines 
in fig. 8, but not in the adjacent spectral regions. The short (long) wavelength interval 
corresponds to a black (white) hole configuration. The presence of optical horizons leads to 
an enhancement of the emission on the background of the horizonless emission that seems 
to decay from the ultraviolet to the infrared. Modes moR and loL exhibit clear horizon 
emission profiles between the black and orange dashed lines and their intervals of emssion 
are indicated by arrows. Over the visible range, emission from moR dominates. Fig. 8 
also shows traces for lower nonlinearities. As expected, the spectral density decreases and 
the invervals of optical horizons, associated with strong emission, become smaller. The red 
line corresponds to zero moving frame frequency u'] no major spectral features seem to be 
associated with this position. 

The quantum state at the output is expected to be a two-mode squeezed vacuum state, 
if two modes only were involved. However, for each moving frame frequency, each mode can 
couple to up to five positive norm and three negative norm modes. Thus, we expect the hnal 
quantum state on the optical branch to be in a partly mixed state across the optical modes. 
Yet, coupling between particular mode pairs seems to dominate in parts of the spectrum, in 
particular within the optical branch. Further characterization of the exact state emerging 
is needed. 


VI. CONCLUSION 

Emission of light from the quantum vacuum has recently attracted attention [1, 14, 19, 32] 
due to the high emission rates expected in the presence of an analogue horizon and the hope 
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to observe analogue Hawking radiation. Using a fully analytical model, introduced in [1], 
the spontaneous emission of light from a moving dielectric boundary, a refractive index 
front (RIF), can be calculated. The model can be used to study any dispersive medium 
phenomenologically described by Sellmeier equations. 


In this work we have focussed on the case of SLIs on both sides of the RIF. We calculate 
the emission for a system based on bulk fused silica. Our approach is interesting, because 
we can show how the emission is structured into intervals with black hole-, white hole, 
and no horizon and how these horizons lead to an increase in the spectral density in the 
laboratory over a general decay. To our understanding, the enhancement of emission due 
to horizons is a common feature of realistic refractive index changes. For index changes 
beyond the damage threshold no enhancement was observed [19]. Additionally, because we 
consider contributions from all modes, we have found that the emission is dominated and 
characterized by a peak in the (laboratory) UV that corresponds to photons emitted in the 
negative norm mode. We show that the emission spectrum from the analogue black hole has 
a unique shape, which is independent of the index change. It is not of the thermal shape 
that is expected from a nondispersive analogue black hole horizon. We also calculate the 
emitted flux and, as a result, the flux increases as a power law of exponent 2.5 for lower 
index steps, as more and more frequency modes appear to participate in the emission at 
the analogue black hole horizon. Whereas previous calculations have mainly focused on the 
Hawking-like emission in the moving frame, these results reveal signatures of the quantum 
vacuum emission important for the experimental observation, such as the nonlinear scaling, 
spectral shape, or mode conhguration. 


This framework forms the basis to future calculations of the revealing quantum correla¬ 
tions between different outgoing modes, correlations of Hawking partners as well as hori¬ 
zonless emission into other modes. This analysis would allow for a better prediction of the 
exact quantum state produced, a key connection to the Hawking effect. The analysis can 
also be extended to more general variations of the refractive index as well as broader choices 
of media, such as waveguides. 
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FIG. 1: Sellmeier dispersion relation, Eq. (2), with three resonances in the laboratory 
frame. There are eight branches (black lines). A contour of ca' is shown in blue. Their 
intersection points indicate the modes of propagation in the medium (red). 
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FIG. 2: Sketch of the refractive index front (RIF) in the moving frame: there are two 
homogeneous regions of uniform refractive index on the left and right of a dielectric 
boundary of height hn. For suitable moving frame frequencies, light on the left (right) of 
the RIF may not (may) propagate away to the right as the step moves at superluminal 
(subluminal) speeds. Light on the left is trapped behind an analogue horizon at C = 0. 
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FIG. 3: Sellmeier dispersion relation of fused silica in a frame moving at a velocity 
u = 0.66 c. Part of the optical branch is shown: branches with positive (negative) 
laboratory frequencies are represented by thick (thin) curves. A curve for zero refractive 
index change 5n is shown in black and for a high change, 6n = 0.02, in orange. Frequency 
intervals corresponding to black and white hole analogue horizons are shaded in orange 

and blue, respectively. 
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FIG. 4: Diagrammatic explanation of the possible mode conhgurations for 
positive-frequency optical modes for various comoving frequencies. Gonfigurations change 
if two of the modes share their group velocity with the RIF (green lines), which occurs at 

the frequencies indicated above and in fig. 3. 
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FIG. 5: Spectrum of emission into the uniquely right propagating mode moR on the right 
side of the RIF. Emission is calculated in the moving frame of velocity u = 0.66 c. The 
number of particles per time and bandwidth, the flux density, is displayed for increasing 
values of 5n (1 * 10“^, 2 * 10“^, 4 * 10“^, 1 * 10“^, 2 * 10“^ (orange dot dashed), 3 * 10“^, 
4 * 10“^, 5.6 * 10“^). The spectrum is confined to a finite interval over which the mode 

moR exists. 
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FIG. 6: Spectrum of emission as in fig. 5. To compare the shapes of the traces, spectral 
densities are scaled such that all traces line up with the 6n = 0.02 trace (orange). Values 

of 6n range from 4 * 10“^ to 5.2 * 10“^. 
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FIG. 7: Emission spectra of each optical mode in the comoving frame 
{u = 0.66 c, 5?7. = 0.02). Emission into mode noL\ purple; loL\ blue dashed; uoL\ green 

dotted; moR\ orange dot dashed. 


21 













10 ^° 

>, 10 ^® 

V) 

C 

cu 

^ I0i6 

03 

l_ 

4 ^ 

U 

cu 

Q. iol4 

13 5 7 

Wavelength A (gm) 

FIG. 8: Emission spectral density in the laboratory frame. At each wavelength the total 
spectral flux density, the number of photons emitted per unit time and unit bandwidth, is 
the sum of contributions from all modes. Emission is concentrated in the UV in a narrow 
spectral peak generated from mode noR. Outside the peak the emission decreases except 
for specral intervals corresponding to black- and white- hole analogue horizons. Spectra are 
calculated for wavelengths above the violet line, beyond which there are no contributions 
from the uppermost dispersion branch. The red line corresponds to ca' = 0 (phase velocity 
horizon). The black and orange dashed lines indicate the interval of the black (white) hole 
mode conhguration for the moR {loL) mode at short (long) wavelengths. 





FIG. 9: Total emission over the subluminal (black-hole) interval into the unique right 
moving mode moR. Estimated photon number for different velocities u (left panel) and size 
of the interval in the moving frame (right panel) as a function of index change Sn. 
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